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ABSTRACT 

As the first step to determine disk mass-to-light ratios for normal spiral galaxies, we present 
the results of an imaging survey for planetary nebulae (PNe) in six nearby, face-on systems: 
IC 342, M74 (NGC 628), M83 (NGC 5236), M94 (NGC 4736), NGC 5068, and NGC 6946. Using 
Blanco/Mosaic II and WIYN/OPTIC, we identify 165, 153, 241, 150, 19, and 71 PN candidates, 
respectively, and use the Planetary Nebula Luminosity Function (PNLF) to obtain distances. 
For M74 and NGC 5068, our distances of 8.6 ± 0.3 Mpc and 5.4±°;4 Mpc are the first reliable 
estimates to these objects; for IC 342 (3.5 ±0.3 Mpc), M83 (4.8 ±0.1 Mpc), M94 (4.4t°;* Mpc), 
and NGC 6946 (6.1 ± 0.6 Mpc) our values agree well with those in the literature. In the larger 
systems, we find no evidence for any systematic change in the PNLF with galactic position, 
though we do see minor field-to-field variations in the luminosity function. In most cases, these 
changes do not affect the measurement of distance, but in one case the fluctuations result in 
a ~0.2 mag shift in the location of the PNLF cutoff. We discuss the possible causes of these 
small-scale changes, including internal extinction in the host galaxies and age/metallicity changes 
in the underlying stellar population. 

Subject headings: distance scale — galaxies: distances and redshifts — planetary nebulae: general 
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1. INTRODUCTION 

Planetary Nebulae (PNe) are dying low-mass 
(M < 8Mq) stars composed of an expanding rar- 
efied shell o f gas surrounding a soon-to-be white 
dwarf fe.g.. iKwok et all 2003 ). The central stars 
of PNe can be among the brightest stellar objects 
in a galaxy with maximum luminos ities above 
~6000 L e (|Vassiliadis fc Wood! Il99l . Due to 
the physics of photoionization, up to ~10% of 
this flux can be reprocessed into a single emission 
line of doubly ionized oxygen at 5007 A, trans- 
forming the o bjects into ve r y bright monochro- 
matic stars (Jacobv 19891 : Jacobv Ez Ciardullol 
19991: iMarigo et all 12004 ) . Moreover, because all 
stars with main-sequence masses between ~1 Mq 
and ~8 M Q eventually evolve into PNe, these ob- 
jects are plentiful in all stellar populations with 
ages between 10 8 and 10 10 yr. 

The above properties make PNe extremely 
useful probes of (1) the i nternal kinematics of 
spiral (iMerrett et al I [20061: ICiardullo etHl 12004 



Dou glas et al.l 120001). elliptic al (jDouglas et al 



2007;_ De Lorenzi et al.| 2008 a nd int eracting 



(|Durrell et al.l 120031: IPeng et all |2004) 



ies, (2) 
Group 



the chemical evolution 



galax- 
of the Local 



e.g., Dopita et ail Il997 : Richer et al 



19981 iMagrini et a l.||2007l), ( 3) the stellar popu- 



popu- 
2006; 



lations of early- type systems (jBuzzoni et al. 
Ciardullol |2006J), (4) the extragalactic distanc e 



scale (e.g.. ICiardullo et alJliobl : ICiardulloll2005l ). 
and (5) the dynamical evolution of galaxy clus- 
ters (iFeldmeier et alj|2005 lArnaboldi et al]|2004 ; 
Gerhard et al.ll2007h . We seek to use PNe as kine- 
matic test particles in face-on spiral galaxies in 
order to use their motions to measure the distri- 
bution of mass in spiral disks. (This analysis will 
be the subject of the second paper in this series, 
hereafter referred to as Paper II.) The first step in 
doing this is to identify large samples of PNe for 
analysis. With these data, we can also measure 
the PNe's photometric properties, use their lumi- 
nosity function to estimate distance, and explore 
how the PNLF behaves in a range of star-forming 
environments. 

In this paper, we describe an [O III] A5007 and 
Ha survey for planetary nebulae in six nearby, 
face-on spiral galaxies. In §2 we present the de- 
tails of our observing runs and image reduction. 
In §3 we explain how we identified our PN candi- 



dates and measured their magnitudes. We discuss 
the discrimination of possible contaminants in §4, 
including how we used photometry and follow-up 
spectroscopy to eliminate H II regions, supernova 
remnants, possible background galaxies, and four 
asteroids from the sample. In §5, we show how we 
determined the completeness limits of our samples 
and also compare our narrow-band selected PNe 
sample in M94 to a sample of planetaries found 
via slitless spectroscopy. In general, we find good 
agreement between the two data sets, although the 
additional information provided by our Ha pho- 
tometry does allow us to eliminate a few compact 
H II regions from consideration. In §6, we derive 
the PNLF distance to each of our six galaxies and 
compare these new measurements to estimates in 
the literature. In all cases, our values are consis- 
tent with previous distance determinations, and 
in three cases, our measurements are significantly 
more precise than those of previous studies. In 
§7 we study the shape of the PNLF, and test for 
differences between various PN subsamples. We 
show that small variations in the shape of the 
PNLF are not uncommon, although these changes 
rarely affect our ability to measure distance. We 
conclude by discussing the possible implications of 
these fluctuations. 

2. OBSERVATIONS AND IMAGE RE- 
DUCTION 

We chose for our study six nearby spiral 
galaxies, all with low recessional velocities (v < 
700 km s~ x ), all relatively face-on (i < 35°), and, 
with one exception, all with large optical angular 
diameters (r > 10'). These galaxies are listed in 
Table [TJ along with their photometric properties 
and metallicities. The two southernmost objects 
in the sample, M83 and NGC 5068, were imaged 
with the Mosaic II camera on the CTIO 4 m 
Blanco telescope. This instrument is composed of 
8 2048 x 4096 SITe CCDs, giving it a pixel scale 
of 0'.'26 pixel^ 1 and a field-of-view large enough 
to encompass the entire galaxy (i.e., 36' x 36' 
or ^10 disk scale lengths for M83 and more for 
NGC 5068). The remaining galaxies were im- 
aged with OPTIC, an orthogonal transfer CCD 
camer a on the WIYN 3.5 m telescope a t Kitt 
Peak (jHowell et al.ll2003t iToiirv et al.ll2002l) . OP- 
TIC, which consists of two Lincoln Lab CCID28 
2048 x 4096 chips separated by 14", produces a 
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pixel scale of 0'.'14 pixel -1 and a field-of-view, 
9.'56 x 9'56, that is slightly smaller than the typical 
size of our galaxies. As a result, our surveys with 
this instrument extended over ~1.2, ~1.8, ~3.5, 
and -7 disk scale lengths for IC 342, NGC 6946, 
M74, and M94, respectively. Although OPTIC is 
designed to track rapid image motion by real-time 
clocking of charge in both the vertical and hori- 
zontal direction, the lack of suitably bright field 
stars prevented us from using this feature. Tabled 
presents a log of our observations. 

For our survey, we followed the standard ob- 
serving procedure for imaging extragalactic plane- 
tary nebulae (e.g., Jacobv et al.ll992t Feldmeier et al 
1997). First we chose a narrow-band filter to iso- 



late the [O III] A5007 emission line at the redshift 
of the target galaxy. For M83 and NGC 5068, this 
filter was custom made, with a central wavelength 
of A c = 5001 and a full-width-half- maximum 
(FWHM) of 46 A in the converging f/2.87 beam 
of the telescope. For the remaining galaxies, stan- 
dard Kitt Peak ~50 A FWHM filters were used. 
These filters are listed in Table [2j models for their 
transmission curves at dome temperature are dis- 
played in Figure [1] along with the wavelength of 
[O III] A5007 shifted to the systemic velocity of 
each target galaxy. 

Each galaxy was imaged via a series of 30 minute 
exposures, with dithering sequences sometimes 
used to fill in the gaps between chips. For M83 
and NGC 5068, corresponding 1 minute exposures 
through a standard V filter served as off-band 
images; for the WIYN observations, continuum 
observations were obtained via a set of 7.5 min 
exposures through the WIYN 32 filter (A c = 5288, 
FWHM = 230 A). In general, these off-band ex- 
posures were scaled to go at least 0.2 mag deeper 
than their on-band counterparts to avoid false 
detections at the frame limit. To help discrim- 
inate PNe from (lower-excitation) H II regions, 
each galaxy was also imaged in Ha. For M83 and 
NGC 5068, the Ha filter had a central wavelength 
of A c = 6563 and a FWHM of 80 A, while for the 
other four galaxies, A c = 6555 and FWHM = 46 A. 
Finally, for IC 342, M74, M83, and NGC 5068, im- 
ages were also taken through a standard R filter. 
These served as off-band images for the Ha expo- 
sures. 

All i mages were r educed using packages in 
IRAF 1 (|Valded Il998h . The Mosaic II expo- 
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Fig. 1. — Transmission curves for the [O III] A5007 
filters used in this survey. These curves repre- 
sent the filter at the outside temperature and in 
the converging beam of the telescope. The wave- 
lengths of [O III] A5007 at the systemic velocity of 
each galaxy are marked. 



sures of M83 and NGC 5068 were reduced en- 
tirely within the mscred package: the data were 
trimmed, bias subtracted, and flatfielded via 
cedproc, placed on an astrometric system by ap- 
plying mscmat ch to ~2800 stars from the USNO- 
B1.0 catalog ([Monet et al. l2003h . projected onto 
the tangent plane with msc image, sky subtracted 
with msc sky sub, aligned with mscimatch, and co- 
added (with cosmic-ray removal) with mscstack. 
During this process, special care was taken to de- 
fine the images' world coordinate system, as we 
found it necessary to run msccmatch separately on 
each individual readout section of each chip (i.e., 
16 runs per exposure) to obtain a proper stack- 
ing and an astrometric precision of < 0'.'3. We 
found the NOAO Deep Wide-Field Survey MO- 
SAIC Data Reductions website 2 to be very useful. 

The OPTIC images of IC 342, M74, M94, and 
NGC 6946 were reduced using a combination of 
IRAF's imred and mscred packages. First, all 



X IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the National Science Foundation. 

2 http: / /www. noao.edu/noao/noaodeep/ReductionOpt /frames. html 
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individual exposures (including biases and flats) 
were converted to Multi-Extension FITS (MEF) 
format using mkmsc; this enabled us to use the 
mscred package to correct for the four different 
overscan regions and the 14" gap between chips. 
The OPTIC reduction notes 3 available off the 
WIYN website were very useful for this first step. 
Next ccdproc was used to trim and flatfield as 
above. Again, special care was needed to de- 
fine the frames' world coordinate system, as each 
readout section (4 per frame) required a man- 
ual initialization with a centrally located star (via 
ecsetwes) and then a refinement (with ccfind 
and ccmap). For M 94, only -20 USNO-B1.0 stars 
( Monet et al.l 2003 ) per quadrant were available 
for the astrometric solution; for M74, this number 
was closer to 40. (Since IC 342 and NGC 6946 
are both close to the Galactic plane, more than 
100 stars per quadrant were available for their 
solutions.) In all cases, our final astrometry at- 
tained a precision of < 0'.'3. Finally, msc image 
was run to project the frames onto the tangent 
plane, linmatch was used to determine the scal- 
ing between exposures, and imcombine was used 
to remove cosmic rays and stack the images. 

3. PLANETARY NEBULA CANDIDATE 
IDENTIFICATION AND PHOTOME- 
TRY 

PN candidates were identified by examining the 
stacked [O III] images, their Ha and off-band 
counterparts, and difference images created by 
scaling, then subtracting the off-band images from 
their on-band frames. PN candidates had to (1) 
have an image profile consistent with that of a 
point source, (2) be present in [O III] but invisible 
in the off-band image (or extremely weak in V), 
and (3) be invisible or weak in Ha (and R when 
available). This latter condition served to elim- 
inate H II regions and supernova remnants from 
the sample (see below). 

Photometric measurements of all the candi- 
dates in [O III] A5007 and Ha were pe rformed 
with the IRAF version of DAOPHOT (jStetson 
1987t IStetson et al.lll990h . First, the point spread 
function (PSF) fitting routines in allstar were 
used to define the PN magnitudes on the differ- 
ence images relative to those of field stars on the 



on-band frames. These results were then checked 
by using DAOPHOT's substar option to subtract 
off a scaled-PSF representation of each PN from 
its position on the on-band and difference frames. 
Although the majority of these subtractions pro- 
duced little or no residual, a few objects, usually in 
or near the galaxy's spiral arms, were found to be 
significantly under- or over-subtracted, sometimes 
by —0.5 mag. The cause of this error was an in- 
correct estimate of the background: in regions of 
active star formation, bright emission from H II 
regions, supernova remnants, and the diffuse ISM 
can play havoc with the underlying "sky," partic- 
ularly in Ha. When this occurred, the PN mag- 
nitudes were manually adjusted until the residu- 
als on the star-subtracted images appeared rea- 
sonable. 

Next our relative [O III] A5007 and Ha mag- 
nitudes were placed on the standard AB sys- 
tem by comparing large-aperture measurements 
of our field stars to similar measurements of spec- 
trophotomctric standards. For our southern hemi- 
sphere measurem e nts, th ese standards came from 
Stone k Baldwin] (|l983h : the stars in the north 



were drawn principally from the list of I Stone 
(|1977l ). Finally, the narrow-band AB magni- 
tudes were converted to monochromatic fluxes 
using models for the filter transmission curves 
in the converging beams of the telescopes (Fig- 
ure []}, knowledge of the galaxies' systemic ve- 
locities (listed in Table [l}, and the photometric 
procedures fo r emis sion-line objects described by 
IJacobv et ail (|l987t) . In the case of our [O III] 
measurements, these fluxes were converted to 
[O III] A5007 magnitudes using 



3 http: / /www. wiyn.org/OPTICreduce.pdf 



™ 5 007 = -2.5 logF 50 o7 - 13.74, (1) 

wher e F5007 is given in ergs cm~ 2 s _1 ( Jacobvl 
1989). At this time, we also used the photometric 
uncertainties derived by DAOPHOT to determine 
the mean photometric error of our measurements 
as a function of [O III] A5007 brightness. This pho- 
tometric error function, which is given in Table [3j 
was used to help determine the PNLF distances of 
§4- 

4. DISCRIMINATING PNE FROM CON- 
TAMINANTS 

Planetary nebulae are not the only sources that 
are detectable via our on-band/off-band survey 
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technique. Deep images through an [0 III] A5007 
filter will also detect H II regions, supernova rem- 
nants, background [O II] emitting galaxies (at 
z ~ 0.34), high-redshift Lya sources (at z ~ 3.12), 
and even asteroids. We discuss each of these be- 
low. 

4.1. H II Regions 

At the typical distance of our galaxies (~5 Mpc), 
1" corresponds to ~25 pc. Consequently, at 
groundbased resolution, compact H II regions will 
be unresolved, and can be misclassified as plane- 
taries. Fortunately, a quantitative way exists to 
eliminate most of these contami nants using two 



of a P N's brightest emission lines. iCiardullo et al 



(|2002j ) have shown that bright PNe populate a dis- 
tinctive region in [O III]-Ha emission-line space. 
Specifically, when R = 7(A5007)//(Ha+[N II]) is 
plotted against absolute [O III] magnitude, true 
PNe occupy a cone defined by 



4 > logi? > —0.37 Mi 



5007 



1.16. 



(2) 



This equation, which is valid over the top ^4 mag 
of the PN luminosity function, is equivalent to the 
statement that the ratio of [O III] A5007 to Ha 
depends on the luminosity of the planetary, via 



4 > R > 3.14 



0.92 



(3) 



where L* = 2.4 x 10 ergs s is the absolute 
[O III] A5007 luminosity of a PN at the bright- 
end of the PNLF. Practically speaking, this means 
that PNe in the top ~1.5 mag of the [O III] lumi- 
nosity function always have [O III] A5007 brighter 
than Ha. This contrasts with the line ratios of 



Ha much brighter than [O III] (e.g., 


Shaver et al. 


1983: Kniazev et al. 2005: Peha et al. 


2007fk 



For three of our galaxies, our narrow-band pho- 
tometry was sufficient to locate each PN candi- 
date in the [O III] A5007-Ha emission line dia- 
gram. Candidates which failed to satisfy the crite- 
rion given by equation ([2|) within the uncertainties 
were excluded from further consideration. Unfor- 
tunately, due to clouds and poor seeing, the Ha 
images of the other three galaxies (IC 342, M74, 
and NGC 6946) did not go deeply enough to de- 
tect the Ha emission from all the possible PNe. 
Specifically, in IC 342 only ~55% of our [O III] 



A5007 point sources were detected in Ha; in M74, 
the Ha recovery fraction was just ~20%, and in 
NGC 6946, none of our bright PN candidates were 
detected in Ha. In the case of IC 342. its near- 
ness (3.5 ± 0.3 Mpc; see §4) and extreme extinc- 
tion (Ay — 1.94 mag) helped to enhance Ha with 
respect to [O III] A5007, somewhat lessening the 
shallow Ha imaging problem. 

For IC 342 and M74, we compensated for the 
lack of Ha photometry using spectroscopy. Ap- 
proximately 68% of the PN candidates in IC 342 
have spectra obtained from the Hydra spectro- 
graph on the WIYN telescope. (A full description 
of these data is presented in Paper II.) For M74, 
the fraction of PNe with spectra is ^65%. Al- 
though these data were obtained for the purpose 
of obtaining radial velocities and only covered the 
wavelength range ~4500 A to ^5500 A, it is pos- 
sible to use them to estimate the Ha flux of each 
candidate. 

To do this, we began by subtracting the con- 
tinuum from each spectrum, and measuring the 
relative line fluxes of [O III] A5007, [O III] A4959, 
and H/3. Appropriately, the ratio of the two oxy- 
gen lines was (on average), three-to-one, indicating 
that the response function within this small wave- 
length range was approximately flat. We then ex- 
amined the throughput of the spectrograph using 
observations of a standard star. Again, the data 
indicated only a slight (~8%) decrease in efficiency 
between the wavelengths of [O III] A5007 and H/3. 
(This is close to the ~5% drop expected from the 
blaze function of the grating.) Using this knowl- 
edge, we derived the true [O III] A5007-H/3 ratio 
for each object, then scaled these values to Ha, 
using an estimate of fore gr ound Galactic e xtinc- 
tion (jSchleeel et al.lll998l) . 



Cardelli et all (|l989h 



reddening law with Ay = 3.1, and an expected in - 
trinsic Ha to H/3 ratio of 2.86 (|Brocklehursdll97l[ ). 
For cases where the Ha flux was not detected pho- 
tometrically, this spectroscopic measurement was 
used to estimate R and eliminate compact H II 
regions from the sample. We note that this tech- 
nique may overestimate the [O III]-to-Ha flux ra- 
tio for PNe that are heavily reddened by circum- 
stellar dust. However, these objects are likely to 
fall well below the PNLF's bright-end cutoff, and 
not be an important part of the sample. 

In the extreme case of NGC 6946, where no 
spectroscopy was available, the best we could do 
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is place an upper limit on the Ha brightness of 
our PN candidates. To do this, we began by using 
the brightest point source magnitude returned by 
allstar during our analysis of the Ha frame. We 
then scaled the frame's PSF to this magnitude, 
and used substar to subtract off this "bright- 
est star" from the location of each PN. We then 
visually inspected the resulting image. As ex- 
pected, due to differences in the underlying back- 
ground, this procedure occasionally resulted in an 
oversubtraction. When this occurred, we manu- 
ally adjusted the limiting allstar magnitude un- 
til the subtraction residual showed no sign of a flux 
deficit. This magnitude was then used to define 
the upper limit to the object's Ha flux. Unfortu- 
nately, even with these limits, it was impossible 
to discriminate PNe from H II regions with sizes 
under —25 pc. Thus, all the PN identifications in 
this galaxy must be considered tentative. 

4.2. Supernova Remnants 

Although compact supernova remnants (SNRs) 
are much rarer than PNe, they may present a 
problem at the very bright end of the PNLF. This 
is particularly true for the more distant galax- 
ies, where groundbased ~1" seeing limits our spa- 
tial resolution to —25 pc. The difficulty is well- 
illustrated by iRoth et aL ( 20041 ). who used spa- 
tially resolved spectrophotometry to show that ob- 
jec t 276 in the list of M 31 bulge PNe presented 
by ICiardullo et al" ( 1989f ) is actually a supernova 
remnant. Fortunately, most SNRs have line ratios 
that clearly distinguish themselves from planetary 
nebulae. For example, of the —60 M31 super- 
nova remnants surveyed by iGalarza et al.l ( 19991 ). 
all but one has [O III] A5007 fainter than Ha. 

More specifically, we can compare the photo- 
metric properties of PNe directly with those of 
supernova remnants using the optical sa mple of 
M83 SNRs created bv lBlair k Lontd (|2004) . Their 
catalog contains 71 supernova remnants, selected 
on the basis of strong [S II] A6716, 6731 relative 
to Ha. Of these, only three fall within 3'.'5 of one 
of our PN candidates, and in each case, there is 
a better, more distinctive SNR candidate nearby. 
None of the other 68 SNRs are within 10" of a PN. 
iBlair &: Lonel ( 2004 ) also presented 18 additional 
[O III]-bright objects, which they believe might be 
extremely young SNRs. Interestingly, all of these 
sources are significantly fainter than the prototype 



of the class, the Milky Way remnant Cas A. They 
are, however, comparable in brightness to the well- 
known SMC remnant, IE 0102-7219, and in the- 
ory, such obje cts can be con f used w ith PNe. Sev- 
enteen of the IBlair & Long] ||2004 ) [O III]-bright 
sources are more than 5" from any of our candi- 
dates and thus not on our list. The remaining 
source (08) is coincident with our M83 PN candi- 
date 10, which satisfies our selection criteria in ev- 
ery respect. Though we targeted this object with 
the Hydra spectrograph, the resulting spectrum 
has very strong Ha emission. While there is vir- 
tually no Ha emission at the location of PN 10, 
there is a strong Ha emitter within 2". We be- 
lieve our spectrum is contaminated by this nearby 
Ha emitter due to a slight misplacement of the 
fiber by Hydra. 

In addition, there are 17 historical supernovae 
in our target galaxies 4 . An examination of our 
images demonstrates that none of the remnants 
from these sources are in our PN sample. We 
did find an emission line object within l'/5 of 
the location of both SNe in M74 (SN 2003gd and 
2002ap) but found nothing of interest at the lo- 
cations of SNe 1983N, 1968L, and 1957D in M83 
and only well-resolved emission-line sources within 
5" of the locations of SNe 1950B, 1945B, and 
1923A (also in M83). Of the nine historical SNe 
in NGC 6946, we found emission line objects near 
three (1980K, 1948B, and 1939C), some H II re- 
gions and blue stars near three others (2002hh, 
1969P, and 1917A), and only star-like objects 
around the remaining three (2008S, 2004et, and 
1968D). 

4.3. Background Galaxies 

Star-forming galaxies at z ~ 0.34 have their 
[O II] AA3726,3729 doublet redshifted into the 
bandpass of our [O III] A5007 filter. Consequently, 
it is theoretically possible for small (< 4 kpc) star- 
bursting objects to be mistaken for planetaries. 
In practice, however, such a misidentification is 
extremely unlikely. At z ~ 0.34, our survey vol- 
ume is -8500 Mpc 3 for the CTIO 4 m Mosaic II 
data, and only —700 Mpc 3 for the galaxies ob- 



4 The Central Bureau for Astronomical Tele- 
grams (CBAT) List of Supernovae website 
( http: / /cfa-www. harvard.edu/iau/lists/Supernovae. html ) 
was used to obtain the locations of the 17 historical SNe 
described here. 
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served with OP TIC. Acco r ding; t o the luminos 
i ty fu nction s of | Hogg et aT (1998), iGallego et al 



(|2002h . and lTeplitz et all (12003). this implies that 



our M83 and NGC 5068 fields contain less than 
—200 [O II] galaxies above our completeness limit, 
and that our OPTIC images contain no more than 
~20 [O II] emitters. Less than 2% of these ob- 
jects will have the extreme rest frame equivalent 
widths (>60 A) needed to appear on our on-band 
frame, b ut be completely invisible in the con- 
tinuum ( Hogg et al.l Il998). Moreover, many of 
these objects will be (at least) marginally resolved. 
Thus, the number of background [O II] interlopers 
in our sample must be negligible. 

In principal, Lya emitting galaxies (LAEs) are 
another source of concern. These objects can have 
extremely large (>100 A) observer-frame equiv- 
alent widths, and, in the absence of deep con- 
tinuum images, can easily be mistaken for PNe. 
Large numbers of such objects have been discov- 
ered in deep [O III] A5007 surveys of V irgo's in- 
tracluster space ( Feldmeier et all 2004 ). and one 
has even been found ma squerading as a PN in the 
M51/NGC 5195 system (|Durrell et alJl2003h ■ For- 
tunately, the surveys described here do not reach 
the flux levels needed to sample t his population. 



None of the 162 LAEs found by iGronwall et al 
(|2007h in their deep 0.28 deg 2 survey of the Ex- 



tended Chandra Deep Field South would have 
been detected on our images of IC 342, M94, 
NGC 5068, and NGC 6946, and only their two 
brightest sources (one an AGN) would have been 
seen in M83. Our survey of M74 does go deep 
enough to detect high-redshift Lya emitters, but 
this galaxy was surveyed with the small field-of- 
view of OP TIC. Consequ e ntly, the luminosity 
function of IGronwall et al.l ( 2007 ) predicts that 
only ~1 of our PN candidates should be a LAE. 
The spectroscopy of Paper II confirms this result: 
Lya emitters are not an important source of con- 
tamination in our survey. 

4.4. Asteroids 

Since both M74 and M83 are located at low 
ecliptic latitudes (—5° and ~18°, respectively), 
one might expect to find asteroids in the images 
of these two galaxies. Asteroids show up as trails 
on the long on-band exposures since they move 
quickly (~0'.'4 min -1 ) across the field of view and 
as nearly point sources on the short off-band ex- 



posures. This makes them very easily distinguish- 
able from PNe. We found no asteroids on our 
~10' x 10' M74 images but found four on our 
— 37'4 x 36f 5 images of M83. The asteroid trails 
in the stacked on-band image correspond to as- 
teroid sky position angles and speeds of ~323° at 
~0'/33 min" 1 , -260° at ~0'/46 min" 1 , -312° at 
~0'/30 min- 1 , and -299° at ~0'.'26 min" 1 , from 
brightest to faintest. A cursory check of the IAU 
Minor Planet Center website suggests the bright- 
est asteroid on our images might be minor planet 
(16522) Tell, which moves ~0'/4 min" 1 with a PA 
of 325? 1 and is expected to have been ~45" from 
the observed location on each of our four frames 
and have virtually the same motion. However, a 
more detailed analysis would be necessary to ver- 
ify this possibility. 

5. COMPLETENESS OF THE FINAL 
SAMPLES 

After the application of the procedures of §4, 
we were left with the following numbers of strong 
PN candidates in each galaxy: 165 in IC 342, 153 
in M74, 241 in M83, 150 in M94, 19 in NGC 5068, 
and 71 in NGC 6946. The positions of these ob- 
jects are displayed in Figure [H Figure [3] shows 
the location of the objects in the [O III] A5007- 
Ha+[N II] emission-line space. For the latter di- 
agram, we have plotted the photometric line ra- 
tio whenever possible; in cases of unreliable Ha 
photometry, either the spectroscopic estimate or 
a lower limit has been displayed. A table of PN 
positions, fluxes, line-ratios, and radial velocities 
will appear in Paper II. Note that it is very diffi- 
cult to find PNe near the centers of galaxies and 
within the spiral arms, due to the bright stellar 
background and diffuse emission associated with 
regions of active star formation. Since the ease of 
PN discovery increases as one moves to the lower- 
surface brightness outer regions of the galaxies, 
many of our PN candidates are at large galacto- 
centric radii: for example, the percentages of PNe 
beyond 2 y-band scale lengths are ~43% in M74, 
-93% in M83, and -73% in M94. 

The next step in our analysis was to determine 
our level of photometric completeness. To do this, 
we performed a series of artificial star experiments 
using the addstar task within DAOPHOT. First 
we examined each image and found the locations 
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Fig. 3.— The flux ratio between [O III] A5007 and 
Ha+[N II] for the planetary nebula candidates of 
our six galaxies. These ratios have been corrected 
for foreground Gal actic reddening using the ex- 
tinction estimates of lSchlegel et al.l ( 1998T ). but the 
effects of internal and circumstcllar reddening arc 
still present in the figure. The dashed lines in- 
dicate values determined from [O III] A5007-H/3 
spectroscopy; arrows indicate lower limits. The 
black cone-shaped outline illustrates the locus of 
Local Group PNe (|Ciardullo et al.l l2002h . Since 
the vast majority of H II regions have Ha brighter 
than [O III] A5007, a diagram such as this is ex- 
tremely useful for discriminating PNe from low- 
excitation interlopers. 

where [O III] emission from H II regions, super- 
nova remnants, and shocks made PN detections 
problematic. These regions were masked out and 
excluded from further analysis. At this time, we 
also excluded bad CCD columns, detector edges, 
and the regions between the two CCDs of OPTIC. 
This removed -25% of the disk area of M74, M83, 
and M94 within the galaxies' inner three disk scale 
lengths, and —16% of the disk area of IC 342 and 
NGC 6946 within one scale radius. (For the small 
galaxy NGC 5068, no exclusion regions were nec- 
essary.) Next, we scattered a small number of ar- 
tificial stars of varying [O III] A5007 brightness 
throughout our on-band image. By adding no 
more than —10% of the original number of de- 



tected planetaries to each frame, we ensured that 
i mage crowding was not a factor in our analysis 
(|Fleming et al.ll995h . Then we searched the frame 
for these artificial stars using the daofind task 
within DAOPHOT and a threshold value of three 
times the background sky. Each recovery (and 
non-recovery) was noted, and the process was re- 
peated several thousand times until the complete- 
ness fraction, as a function of [O III] A5007 mag- 
nitude and position in the galaxy, was well deter- 
mined. Finally, we fit this function to an equation 
of the form 



/(m) = 1/2 



f3(m - m Um ) 
y/1 + f3 2 (m- m Hm ) 2 



(4) 



([Fleming et al.l Il995l ) . For purposes of this pa- 



per, we define our limiting magnitude for PN de- 
tections to be the 90% completeness limit. This 
number is listed in Tabled 

5.1. Comparison with Previous PN Sam- 
ples: M94 

One of the galaxies in our sample has previously 
been surveyed for planetary nebulae. M94 was 
observed with a dual-beam, "counter-dispersed 
imaging" spectrograph, designed to create Ha im- 
ages with one arm, and measure PN velocities via 
[O III] A5007 slitless spectroscopy with the other 



( Douglas et al. 2000[) . Although PN photometry 



through an instrument such as this is difficult and 
subject to a range of systematic errors, it is still 
meaningful to compare the two samples. 

Douglas et al. (2000) studied two regions along 



M94's major and minor axes, and found 53 and 14 
PN candidates, respectively. Their fields are dis- 
played in Figure [2] Of the 67 PN candidates found 
via slitless spectroscopy, 44 are located within 5" 
of one of our candidates, 16 are likely H II regions 
and not included in our sample, one is outside 
our field-of-view, one appears to be a duplicate, 
and five are below the limit of our completeness 
and not present in our sample. Conversely, 22 of 
the PNe found in ou r survey are not cataloged by 
Douglas et al. (2000). Of these, eight are near the 
edge of the spectrograph's field, and thus may have 
been missed for that reason. Of the remaining 14, 
two are within —0.5 mag of the bright end of the 
PNLF, nine are between 1 and 1.5 mag down the 
luminosity function, and three are extremely faint. 



S 



A comparison of the objects' [O III] magnitudes 
shows generally good agreement. As described 
above, 44 PNe were detected in both surveys. Of 
these, four have slitless spectroscopy magnitudes 
that are more than 0.8 mag fainter than our mea- 
surements, and two have no measured magnitude 
at all. However, the remaining 38 have magnitudes 
that are consistent with our own. The scatter be- 
tween the two data sets is a — 0.22 mag, and the 
systematic zero-point offset between the two sam- 
ples is only 0.019 ± 0.036 mag, with our values 
being slightly fainter. 

Finally, when we examine the objects that we 
tagged as likely H II regions, we again find good 
agreement. Of the 16 sources, 10 were classified by 
IDouglas et all (|2000l ) as H II regions based on their 
extremely brigh t [O III] magni t udes. Conversely, 



one object that Douglas et al. ( 2000h considered 



too bright to be a PN is classified as a planetary 
in our sa mple. This disagreem ent arises from the 
fact that IDouglas et alj (|2000h assumed a distance 
of 6 Mpc in their analysis. As we describe below, 
our distance to the galaxy is substantially smaller 
than this value. 

Of the 44 candidates in common between the 
two samples, we have spectra of 37. A comparison 
of the velocities from these two samples will be 
included in Paper II. 

6. PNLF DISTANCES AND COMPARI- 
SON TO EXISTING ESTIMATES 

The Planetary Nebula Luminosity Function is 
an excellent standard candle: observations in —50 
elliptical, spiral, and irregular galaxies demon- 
strate that the bright-end of this function is well 
represented by an equation of the form 



N(M) oc e 



0.307M 



1 



3(M*-M) 



(5) 



where M* is the absolu te magnitude of the 



most luminous pla netary (jCiardullo et al. 1989; 



Jacobv et al. 19921 ). An external calibration us 



ing galaxies with well-determined Cepheid dis- 
tances defines the nominal value of M* to be 
— 4.471q q3, although this cutoff appears to fad e 
in low-metallicity objects ( Ciardullo et al. 20021 ). 

To derive the PNLF distances, we began by 
convolving the empirical PNLF of equation ([5|) 
with the photometric error functions defined in 



Table [3] We then fit those PNe brighter than 
our 90% completeness limit to the convolved 
function via the meth od of maximum likelihood 
(ICiardullo et al.lll989h . (Since the PNLF of M83 
appears to deviate from this law at faint mag- 
nitudes, we truncated its fit at 7715007 = 25.5.) 
Foreground Galactic extinction was taken into 
account by adopting the e xtinction estimates o f 
Schlegel et all <|l998h and a ICardelli et all (| 19891 ) 
reddening law, i.e., A 500 7 = 3A7E(B - V). This 
analysis yielded the galaxies' distances, along with 
their formal fitting errors. To these errors, we then 
added (in quadrature) the possible systematic un- 
certainties associated with our photometric stan- 
dard star measurements (—0.02 mag), the filter re- 
sponse function (—0.02 mag), and the foreground 
extinction estimate (0.16 x A5007; see Table [4}. 
Finally, a Kolmogorov-Smirnov test was used to 
determine whether the best-fit empirical function 
was indeed a valid representation of the observed 
data. In all cases, this is confirmed to be true: the 
worst fit is for M74, which has a strange dearth of 
PNe between -0.2 and -0.4 mag below the PNLF 
cutoff (see §7). However, even for this galaxy, the 
null hypothesis can only be excluded at the —90% 
confidence level. Figure 0] displays the observed 
luminosity functions for each galaxy, along with 
the best fit models. 

We note that our analysis of the PNLF ig- 
nores the effects of internal extinction. This 
omission is just i fied f or several reasons. First, 



Giovanelli et al.l (|1994T ) have found that Sc galax- 



ies have very little dust extinction beyond two 
or three scale lengths of the galactic center. As 
we mentioned above, many of our PN candidates 
in M74, M83, and M94 have gal actocentric radii 
larger than this limit. Moreover. iFeldmeier et al 



(|19970 modeled the effects of internal extinction 
on the PNLF by assuming that the vertical scale 
heights of dust and PNe in other spiral galaxies 
are similar to those observed in the Milky Way. 
They found that the larger scale height of the PNe 
make the PNLF method relatively insensitive to 
internal extinction. Specifically, derived PNLF 
distances are never more than —0.1 mag less than 
their nominal value, even when the total amount 
of ^5007 extinction is large. In theory, this ques- 
tion can be addressed directly, via spectroscopy 
of the PNe's Ha and H(3 emission lines, though 
the faintness of H/3 and the effects of atmospheric 



9 



PTJTTTTJTTTTJTTTTJTTT 

IC 342 



10 



I | I I I I | I I I I | I I I I | I I I b 
_ S3 

M74 




26 27 24 25 26 

Apparent X5007 magnitude 

Fig. 4.— The [O III] A5007 PNLFs of the six 
galaxies analyzed in this study. The curves rep- 
resent the best-fitting empirical PNLFs convolved 
with the photometric error function and shifted to 
the most likely distance of the galaxy. No correc- 
tions for reddening have been applied here. The 
open circles indicate points fainter than the com- 
pleteness limit. 



dispersion make the observations challenging. 

In the sections below, we compare our PNLF 
distances to values determined from the anal- 
yses of Cepheids, the Tip of the Red Giant 
Branch (TRGB), Surface Brightness Fluctua- 
tions (SBF), the Brightest Blue Super Giants 
(BBSG), the Brightest (Red and Blue) Super- 
giants (BSM), and the Expanding Photospheres 
of Supernovae (EPM). Though some of the galax- 
ies also have Tully-Fisher measurements, the fact 
that we selected our targets to have low inclina- 
tion makes these determinations unreliable. Fig- 
ure [5] and Table [5] compare our distance mod- 
uli to measurements in the literature. Note that 
the zero point for all our PNLF distance s is tie d 
to that of Cep heids via ICiardullo et all (|2002l ). 
Thus, following iFreedman et al.l (|200ll ). our dis- 
tance scale assumes an LMC distance modulus of 
M) = 18.50. 



im 



6.1. IC 342 



IC 342, properly dubbed the Hidden Galaxy, 
would be one of the brightest galaxies in the sky 
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Fig. 5. — Final PNLF distances compared to ex- 
isting distances in the literature. In most cases, 
only the distances with the smallest uncertainties 
are shown, but a few measurements with large un- 
certainties and one without quoted errors are pro- 
vided for comparison. See Table 5 for more infor- 
mation on the distances, including the reference 
for each value. 



except for its location near the Galactic plane 
(b 11 ~ 11°) and the large amount of foreground 
extinction (E(B - V) = 0.558). The galaxy 
is part of the Maffei group, a system that gar- 
nered some attention due to the suggestion of 
a past gravitational encounter with the Local 
Group ( Zheng et all Il99lt IValtonen et"al] Il993l: 



Peebles 



1994). More recently, multiple studies 



have presented evidence against this possibilit; 
dKrismer et al. 1995 : iKarachentsev et "aTl 12003 



Fingerhut et al.l 120071 ) . Although the group of 



~16 galaxies was originally thought to have con- 
tained only two large galaxies (IC 342 and the 
giant elliptical Maffei 1), a recent analysis by 
iFingerhut et al. ( 2007h has resulted in the addi- 
tion of another member, the heavily extinguished 
(Ay ~5.6 mag) spiral Maffei 2. 

A number of different distance indicators have 



been used for IC 342. ISaha et al.l (|2002h used two- 
passband observations of 20 Cepheids to estimate 
a distance modulus of (to— M)o = 27.58±0.18 and 
a total reddening of Ay — 2.01. With an updated 
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period-luminosity relation and metallicity mea- 
surement, this numb er increases slightly to (to — 
M) = 27.64 ±0.12 (jFingerhut et al.ll2007t ). For 
comparison, the magnitudes of the three brightest 
blue and red supergian ts gives a distance modulus 
of (to - M) ~ 27.29 dKarachentsey fc Tikhonov 
19931 revised bv iKarachentsev et alJ (|l997t )). 

Out of the 165 PNe found in IC 342, the 65 in 
the top —0.9 mag of the PNLF qualify as our com- 
plete sample. Our application of the maximum 
likelihood method then gives a true distance mod- 
ulus of (to - M) = 27.72 ± 0.17 {3.50±oit M P C )- 
Note that the error in our distance is completely 
dominated by the uncertainty in the foreground 
extinction. (Without this term, our error would 
on ly be —0.05 mag ! ) Sinc e the extinction derived 
bv lFingerhut et all tj2007h is completely consistent 
with that determine d from the COBE satellite 
( Schlegel et al. 19981 ) . we adopt their —0.16 mag 
uncertainty in our calculation as well. 

Figure [5] compares our distance estimate to that 
of the Cepheids and supergiant stars. Our value 
is statistically identical to the former value, but 
substantially larger than that associated with the 
galaxy's brightest stars. 

6.2. M74 

The Perfect Spiral, M74 (NGC 628), is per- 
haps the most face-o n galaxy in the nearby un i- 
verse, with i ~ 6.5° (jKamphuis &: Briggsl ll992). 
This grand-design system dominates a small 
group which includes th ree dwarf irregular ob- 



jects ( Sharina et al.fl99^ 1. and play s host to both 



our n earest Type Ic "hypernova" (Ma zzali et al 



2002) and a n extremely variab le ultraluminous 
x-ray source (jKrauss et al. Idooi). 

Though M74 is well-studied, it lacks a reliable 
distance estimate. Va lues in the literature r ange 
from -7 to -20 Mpc (|Sohn fc Davidgejll996l ) and 
this uncertainty propagates into all aspects of its 
analysis, including luminosity estimates for its su- 
pernovae and x-ray sources. M74's three least 
uncertain distances come from the photometry of 
its three brig htest blue supergian ts ((to — M) = 
29.32 ± 0.40; ISharina et atlll996l ). the photome- 
try of its brighte st supergiants of all colors ((to — 
M) = 29.5 + g |: Irfendrv et al.l[2005h . and the ex- 
pansion velocity-absolute magnitude relation of 
the Type II-P supernova SN 2003gd ((to - M) = 



29.9±g;|; iHendrv et a0l2005h . 

Our complete sample of planetaries consists of 
82 objects out of a total population of 153. By 
fitting the top —1.1 mag of the luminosity func- 
tion, we obtain a distance modulus of 29.671q'q7, 
or 8.59+Q27 Mpc. Though this distance differs 
by only —10% from the values obtained from the 
Type II-P supernova and the galaxy's supergiants, 
it has considerably less uncertainty. 

6.3. M83 and NGC 5068 

M83 (NGC 5236), the barred Southern Pin- 
wheel, is one of the two central galaxies of the 
nearby Cen A (NGC 5128) group, an —87 galaxy 
complex which also includes NGC 5068. Though 
often viewed as a single system. IKarachentsev et al 
(|2002j, 120071 ) have recently shown that the group 
actually contains two spatially separated associ- 
ations, one centered on the interacting elliptical 
NGC 5128, and the other clustered about M83. 
M83 itself is one of the great supernova factories 
in the sky, with six recorded events in the last 
century. 

To date, the two most accurate distances to 
M83 have come from Cepheids and stars at the tip 
of the red giant branch. The former method, based 
on V and J-band observations of 12 Cepheids with 
the VLT, giv es a distance modulus of (to — M) n = 
28.32 ± 0.13 dThim et all 120031 : ISaha et allEloi k 
the latter, which comes from HST F606W and 
F814W photometry of halo stars, yields (to — 
M) = 28.56 ± 0.18 (jKarachentsev et al.|[2007h . A 
group distance, determined via TRGB measure- 
ments to 10 of the galaxy's surrounding dwarfs, is 
(to - M) = 28.40±°;°5- 

Our PN survey with the CTIO 4-m telescope 
detected 241 PNe in M83; 164 of these fall in 
the top —1.3 mag of the PNLF and comprise the 
complete sample before the luminosity function 
turns over. Our maximum-likelihood fit to the 
data yields a distance modulus of (to — M) = 
28.43±^e, or 4.85^3 Mpc. Interestingly, this 
value lies between the Cepheid and TRGB dis- 
tances to the galaxy, and is closer to the mean 
group distance than either of the two direct mea- 
surements. 

For NGC 5068, 11 of our 19 PNe are above our 
completeness limit and in the statistical sample. 
Formally, our fit to these data produces a distance 
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modulus of (to - M) = 28.85±o°g, a value that 
places the galaxy in back of the main Cen A /M83 
complex. However, ICiardullo et al ] (120021 ) have 
noted that for systems more metal-poor than the 
LMC, the PNLF cutoff fades, leading to system- 
atically larger d i stanc es. To correct for this effect, 
Ciardullo et all (|2002h suggest using the t heoreti- 
cal relation derived bv lDopita et al. ( 19921 ) for the 
behavior of M* with metallicity, i.e., 

AM* = 0.928[O/H] 2 + 0.225[O/H] + 0.014, (6) 

where the solar abundance of oxygen is assume d 
to be 12 + log(0/H) = 8.87 (jGrevesse et al.lll996l) . 
With this relation and the observed low oxygen 
abundance of the ga laxy, 12 + log[0/H] ~ 8.32 
(|Pilvugin et al.ll2004h . NGC 5068 moves closer to 
M83, to 5.45to:4° Mpc. This is the first direct 
distance measurement to this object. 



6.4. M94 

One of the many galaxies in the Canes Ve- 
natici I cloud, M94 (NGC 4736) is surprisingly 
isolated with no visible companions down to a B- 
band surface brightness limit of ' 25 m ag arcsec -2 
( Karachentseva fc Karachentsev 19981) within ^3 
degrees or 230 kpc. This Sab spiral has a small, 
but bright bu lge, a Low Ioniz ation Nuclear Emis- 
sion Region (jHeckmanl Il980h . and an extremely 
bright H II ri ng encircling its nucleus at a distance 



of ~1.1 kpc (|Munoz-Tun6n et al.lliool 



The dis tance to M94 is no t well e stablished. 
Although iKarachentsev et "al" ( 2003bf ) obtained 
a value of (to - M) = 28.34 ± 0.29 from HST 
imag es of the T ip of the Red Giant Branch, 



and iTonrv et all (|200lh derived (to - M) 



28.58 ± 0.18 from Surface Brightness Fluctuations 
of the galaxy's bulge, a significantly larger dis- 
tance ((to — M)o = 28.89), based simply on the 
object's Hubble Flow and Hq = 55 km s _1 Mpc -1 
is often used in the literature. 

We identified 150 PNe in M94, with 71 in a com- 
plete sample which extends ~1.3 mag down the 
PNLF. Our distance modulus from fitting these 

Mpc. 



1+0.12 
-0.16 



data is (to - M) = 28.21+£^, or 4.39 
This easily fits within the uncertainty of the 
TRGB measurement but is incon sistent with the 
more distant SBF value. However. lCiardullo et al 



(|2002l ) have shown that the ze ro point of the SBF 
distance scale as described by ITonrv et al. ( 2001 ) 



is systematically larger than that of the PNLF by 
0.30 ± 0.05 mag. (Ha lf of this dis c repan cy has 
since been resolved by IJensen et al. (2003) via a 
scale adjustment to the SBF calibrators.) If this 
correction is taken into account, then the SBF dis- 
tance becomes statistically indistinguishable from 
our value. 

6.5. NGC 6946 

NGC 6946, the Fireworks Galaxy, has a nu- 
clear starburst, significant star formation activ- 
ity throughout its spiral arms, and an impressive 
number of recent supernovae (9 in all) . It is rela- 
tively isolated, having only ~1 low surface bright 
ness objects in its group (IKarachentsev et al 



l2000h . and it is located close to the Galactic plane, 
(b 11 ~ 12°). It therefore falls behind a relatively 
large amount of Ga lactic extinction (Ab ~ 1.475; 
ISchlegel et aTlll998l ). 

There are few rel iable distances to NGC 6946. 
Shar ina et al. I (|l997h used ground-based B and V 



images to identify three of the galaxy's bright- 
est supergiants, and estimated a distance mod- 
ulus of (to - M) n = 28 90 ± 0.18. Similarly, 
Karachentsev et al. (2000) analyzed the bright- 



nesses of blue supergiants in the galaxy's sur- 
rounding dwarfs and inferred (to — M)q = 28.85 ± 
0.15. Both these values im ply a distance t hat is 
smaller than that found bv lSchmidt et al. ( 19921 ) 
by applying the Expanding Photosphere Method 
to the Type II-L supernova 1980K ((to - M) = 
90 29 +0 - 20n ) 

^■^-0.32)- 

Because the spiral arms of NGC 6946 are filled 
with bright H II regions, identifying PNe in this 
galaxy is particularly difficult. As a result, our PN 
sample consists of only 71 candidates, and only 11 
of these are brighter than our 90% completeness 
limit. Moreover, as detailed in §4.1, our Ha ob- 
servations did not extend deep enough to confirm 
the high-excitation nature of these objects, and 
no follow-up spectroscopy was available in lieu of 
this limitation. Still, the observed shape of our 
PNLF is consistent with that expected from the 
empirical model, and our fit of the top ~0.3 mag 
of this function yields (to - M) = 28.92 ± 0.21 
(6.08tp;f6 Mpc). We note that despite the limited 
number of PNe in the fit, our distance uncertainty 
is still dominated by the error associated with the 
galaxy's foreground extinction, ~0.2 mag. We 
also note that, although some of our PN candi- 
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dates may in fact be H II regions, it is likely that 
our PNLF distance will remain largely unaffected. 
Specifically, a careful examination of our images 
suggests that four of our final 71 PN candidates 
might be slightly resolved. If we eliminate these 
objects from our sample, our derived distance to 
NGC 6946 increases by less than 2%. 

7. VARIATIONS IN THE PNLF 

The PNLF distance method has been tested as 
well as any general purpose extragalactic standard 
candle. Despite these tests, a few lingering doubts 
remain. In the HST Key Project 's calibration of 



the P opulation II distance scale, iFerrarese et al 
(|2000h chose not to include the PNLF results 



in their final distance determination, since its 
distances were systematically shorter than those 
found fro m Surface Brigh t ness F luctuations. More 
recently, ISambhus et al. I (|2006h claimed the de- 
tection of two kinematically distinct PN popu- 
lations in the E6 elliptical galaxy NGC 4697, 
with each following its own luminosity function. 
Moreover, surveys of Coma and Virgo have pre- 
sented evidence for the possible existence of over- 
luminous planetary nebulae in intergalactic space 
(|Arnaboldi et al.ll2007l 120081 ). 

The best way to test the robustness of the 
PNLF method is to examine subsets of PNe within 
a single galaxy. Population gra dients in spiral 



and elliptical galaxies are common (jZaritskv. et al 



1994 IPeletier et allll990l) . so if the PNLF is ob- 



served to be constant across the face of a galaxy, 
that presents strong evidence that neither age nor 
metallicity is affecting the luminosity function. 
Unfortunately, to date, only a few PN popula- 
tions have been studie d in this way. Anal yses of 



the PNe across M31 (Mcrrctt et al 



(Magrini et al 



(Ma grini et al 



2000: Ciardullo et al 



2001), M104 dFord et al 



2006), M33 
2004h M81 



1996), 



and NGC 5128 (|Hui et all 119931 ) have found no 
radial change in th e PNLF, while the gradi- 
ents found in M87 dCiardullo et all Il998l ) 



and 



NGC 4526 (jFeldmeier et al.ll2007l) have been at 



tributed to the presence of foreground intracluster 
stars. 

The planetary nebula samples of IC 342, M74, 
M83, and M94 are large enough for us to exam- 
ine the systematics of subsamples. We therefore 
divided the PNe in these galaxies in three dif- 



ferent ways: radially (i.e., by defining an inner 
and an outer sample), by the direction of rotation 
(assuming that all PNe are rotating along with 
the galactic disk), and by the orientation of the 
galaxy (i.e., near side versus far side). The first of 
these divisions should be sensitive to population 
changes within the program galaxy, the second to 
wavelength-dependent variations in the transmis- 
sion of the narrow-band filter, and the third to 
host galaxy internal extinction. Figure [5] illus- 
trates these divisions. 

Figure [6] compares the cumulative luminos- 
ity functions of our various subsamples via a 
Kolmogorov-Smirnov test. In general, there is ex- 
cellent consistency within the galaxies. However, 
as the statistic shows, a few anomalies do exist. 
In IC 342, the PNLF derived from the sample of 
planetaries on the approaching side of the galaxy 
differs from that derived from PNe on the receding 
side at the 90% confidence level. Interestingly, as 
the differential plot of Figure [7] illustrates, both 
sets of objects yield identical values for the PNLF 
distance. However, while the planetaries on the 
receding side obey the luminosity function given 
in equation O the approaching PNe deviate from 
this law with 90% confidence. The PNe of M74 
show a similar effect, except this time the gradi- 
ent lies in the radial direction. Again, the PNLF 
distances inferred from both samples are nearly 
identical, but the detailed shapes of the two func- 
tions differ with 95% confidence. In this case, 
the problem can be traced to a deficit of objects 
~0.3 mag down from to*. Finally, in M83 and 
M94, the PNe on the near side of the galaxy follow 
a slightly different luminosity function than those 
on the far side (at ^90% confidence). In M94, 
this difference does not affect the determination 
of m*, since it is primarily caused by a deficit of 
objects ~0.5 mag below the cutoff in the sample 
of planetaries north of the galaxy's major axis. 
However, in M83, to* is affected: although both 
samples of planetaries are well fit by equation ([5]) , 
the PNLF for objects on the northwest side of the 
galaxy yields an apparent distance modulus that 
is ^0.18 mag farther than that derived from the 
southeast sample. This discrepancy is ~3 times 
the formal error of the fits. 

There are two possible explanations for the in- 
consistencies seen in the figure. The first is inter- 
nal extinction in the host galaxies. In the Milky 
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25 25.5 26 26.5 25.5 26 26.5 27 24 24.5 25 25.5 26 23.5 24 24.5 25 25.5 

Apparent X5007 magnitude 

Fig. 7.— The differential [O III] A5007 lumi- 
nosity function for the four systems where the 
Kolmogorov-Smirnov statistic indicates a change 
in the PNLF across the galaxy. In the case of 
IC 342, M74, and M94, the most-likely distance 
modulus is not affected by these fluctuations. In 
M83 however, the best-fit value for m* changes by 
~0.18 mag. 



Way, the scale height of P Ne is only slightly larger 
than that of the dust dPandev fc Mahral 1 1987 ; 



Ziilstra fc Pottaschlll991l ; iPhillips 20011). Such an 



intermingling can depress the magnitudes of some 
objects, but not others, and create a luminosity 
function very similar to that seen in the inner re- 
gions of M74 or the receding side of IC 342. Of 
course, to explain all the variations seen in Fig- 
ure El one needs the extinction to be patchy, or 
at least have some azimuthal dependence. More- 
over, it is well known that t he amount of dust in a 
galaxy declines wit h radius (jGiovanelli et al 1 ll994t 
iKvlafis et alll200lh . so that if internal extinction 
is truly responsible for these effects, one would ex- 
pect to see the PNLF change more with radius 
than with position angle. 

An alternative explanation involves assuming 
that the PNLF's shape is population dependent, 
and that the differences seen in Figure [H] are due 
to subtle changes in the age or metallicity dis- 
tribution function of the progenitor stars. To 
some extent, such changes should not be surpris- 
ing. We know that the faint end of the PNLF is 
population dependent: while old systems such as 
M31's bulge have exponentially increasing PNLFs, 
star-forming systems such as the Small Magellanic 
Cloud exhibit a "dip" in the lumi nosity function 
start i ng ~2 mag below the cuto ff ( Ciardullo et alj 
2002Hjacobv fc De Marcoll2002f ). In fact, Figure!! 



suggests that the same dip seen in M33 and the 
SMC may be present in M83 as well. However, 
there has never before been any evidence for small- 
scale features that distort at the bright-end of the 
PNLF. 

It is becoming increasingly clear that the term 
"planetary nebula" refers to a heterogeneous col- 
lection of objects. PNe from single stars will have 
luminosities which depend on the m ass of their 
central star ( Vassiliadis fc Woodlll994 ). Since cen- 
tral star mass depends on progenitor mass through 
the initial mass-fina l-mass relation (jWeidemannl 
200rtlWilliarnll2007h . populations of different ages 
will produce different PNLFs. Binary stars also 
produce planetaries via common envelope interac- 
tions, and their contribution to the PN LF may be 
very different from that of single stars (ISokerll2003l 
Moe fc De Marco! l2006h . Finally. Iciardullo et all 
(2005) have proposed that the descendants of blue 
stragglers may contribute greatly to the bright-end 
of the planetary nebula luminosity function. The 
observed shape of a population's PNLF may in- 
volve a complex mixture of contributions from all 
these sources. The resiliency of the PNLF tech- 
nique may simply be due to the processes which 
place a hard upper limit on the emitted flux at 
[O III] A5007. 

8. DISCUSSION 

The results of our analysis demonstrate that 
the PNLF is still an excellent standard candle. 
The distances presented to the six galaxies dis- 
cussed in this paper are in good agreement with 
previous measurements, and the subsample com- 
parisons performed in §7 indicate that there are 
no major population gradients in the PNLF cut- 
off. However, minor unexplained variations in the 
luminosity function do exist, and in some cases, 
this can introduce an additional ^0.2 mag error 
into the distances. At this time, it is impossible 
to determine whether these variations are intrin- 
sic to the stellar population, or caused by inter- 
nal extinction in the host galaxies. The fact that 
these are star-forming galaxies with large amounts 
of dust suggests t hat extinction is the cause, but 
the data from the lSambhus et al" ( 20061 ) survey of 
the elliptical galaxy NGC 4697 and the observa- 



tion of Virgo's intracluster stars (jArnaboldi et al 
20081 ) suggest otherwise. In theory, one can dis- 
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criminate between these (non-exclusive) hypothe- 
ses via spectrophotometry of the objects' Ha/H/3 
ratio. In practice, however, such a study is ex- 
tremely difficult. In a typical bright planetary, H/3 
is ~10 times weaker than [O III] A5007, so obtain- 
ing accurate measurements is a challenge. More- 
over, PNe are subject to circumstellar extinction, 
as well as Galactic extinction. Disentangling the 
two components will require the analysis of large 
samples of objects, in many different galactic en- 
vironments. 
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Fig. 2. — Our [O III] A5007 images of the six galaxies studied in this survey. North is up and east is to the left; the 
positions of the PN candidates are colored such that red, orange, yellow, green, and blue represent objects in the top 
0.5, 1.0, 1.5, 2.0, and 2.5 mag of the PNLF, respectively. The major and minor axes are shown as white lines. For 
IC 342, M74, M83, and M94, the inner ellipse marks the high surface-brightness region of the galaxy within which PN 
identifications were severely incomplete, while the outer ellipse shows the dividing line between our inner an d outer 
samples. For M94, the green rectangles indicate the approximate field of view of the iDouglas et "all (|200Ct ) slit less 
spectroscopy survey, the black crosses indicate candidates in both samples, and the yellow diamonds indicate their 
candidates that are not present in our sample. 



16 



Inner/Outer 



20 
10 




I I I I I I 




- IC 342 




- D m „=0.133 




II I IT I I I 


, I 


25.5 


26 



To ward /Away 
n — i — | — i — i — i — i — 
30 |- IC 342 

D mox =0.313 



Near/Far 



I I I I I I 


1 1. 


- IC 342 




- D mox =0.211 




n iTT i i i 


1 1 " 




24 24.5 25 

Apparent A5007 magnitude 

Fig. 6. — The cumulative [O III] A5007 luminosity functions for subsamples of PNe within four galaxies, with 
the Kolmogorov-Smirnov test statistic marked. In four out of the 12 comparisons, this statistic excludes the 
null hypothesis with more than 90% confidence. The approaching PNe of IC 342, the inner PNe of M74, and 
the southwest PNe in M94 all deviate from the empirical law with more than 90% confidence. Although the 
near-side and far-side samples of PNe in M83 both follow the empirical law, the cutoff magnitude changes 
by ~0.18 mag. 
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Table 3 

Photometric Error versus Magnitude 
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0.155 


0.093 


0.184 






0.180 


26.8 


0.175 


0.110 


0.205 






0.207 


27.0 


0.198 


0.135 


0.229 






0.232 


27.2 


0.225 


0.163 


0.255 






0.269 


27.4 


0.260 


0.192 








0.325 


27.6 


0.306 


0.225 











Table 4 

Possible Systematic Uncertainties (in mags) 





Photometric 


Filter 




Galaxy 


Zero Point 


Response 


Extinction" ' fc 


IC 342 


0.015 


0.010 


0.160 


M74 


0.028 


0.010 


0.038 


M83 


0.015 


0.010 


0.038 


M94 


0.040 


0.030 


0.010 


NGC 5068 


0.010 


0.010 


0.056 


NGC 6946 


0.037 


0.010 


0.191 



References. — a: Schlegel et al. 1998; b: Cardelli et al. 
1989 
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Table 5 

Comparison of Distance Moduli 



Galaxy 



(m - M) D (in Mpc) 



Method 



Reference 



IC 342 
IC 342 
IC 342 



27.64 ± 0.12 
27.72 ± 0.17 
27.29 



3.4 ± 0.2 

3 ™ +0 ' 28 
o.ou„ 26 

2.87 



Revision of WIYN Ccphcids 
PNLF 

Blue Supergiants corrected for H II regions 



HE' 



erhut et al. (2007 



This study 
Karachentsev et al. 



M94 
M94 
M94 



28. 2r 



I +0.06 
-0.08 

28.42 ± 0.18 
28.34 ± 0.29 



4.39~* 



3 + 0.12 
-0.16 

4.8 ± 0.4 
4.7 ±0.6 



PNLF 
SBF 

HST TRGB 



This study 
Tonrv et al. (2001) ; Jensen et al. (2003 ) 



Karachentsev et al. (2003b 



M83 
M83 
M83 
M83 



28.40±° ;04 

28.43ig;SS 
28.32 ± 0.13 
28.56 ± 0.18 



4.8 ± 0.1 
4 85+ ' 11 
4.6 ± 0.3 
5.2 ± 0.4 



HST TRGB of group 
PNLF 

revised VLT Ccphcids 
HST TRGB 



Karachcntscv^^^ ^007) 
This study 
Saha et al. (2006) 
Karachentsev et al. (2007) 



NGC 5068 



28.68" 



0.18 



5.45- 



0.13 



Metallicity corrected PNLF 



This study 



NGC 6946 
NGC 6946 
NGC 6946 
NGC 6946 



28.85 ± 0.15 
28.90 ± 0.18 
28.92 ± 0.21 



2'l.2<) 



+0.20 
-0.32 



5.9 ± 0.4 
6.0 ± 0.5 

k n«+ - 62 

b - U8 -0.56 

7 9+0.7 

'•'-1.0 



Brightest Blue Supergiants of group 

Brightest Supergiants 

PNLF 

Expanding Photospheres Method 



Karachcntscye^^ ^000) 
Sharina et al. (1997) 
This study 
Schmidt et al. (1992) 



M74 
M74 
M74 
M74 



29.67! 



7 +0.06 
-0.07 

29.32 ± 0.40 

c+o.5 



29.5 
29.9 



0.7 
+0.6 
0.8 



8.59"* 



H0.24 
-0.27 
7 Q + l-5 
'' J -1.2 

8.1 ± 2.3 
9.6 ± 2.8 



PNLF 

Brightest Blue Supergiants 
Brightest Supergiants 
SN Il-P expansion rate 



This study 
Sharina et al. (1996) 
Hendry et al. (2005) 
Hendry et al. (2005) 



Note. — Galaxies are listed by increasing distance and distance moduli are listed by increasing uncertainty. 
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